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ABSTRACT 

We present an analysis of the ENEAR sample of peculiar velocities of field and cluster elliptical galaxies, 
obtained with D n -a distances. We use the velocity correlation function, ipi(r), to analyze the statistics 
of the field object's velocities, while the analysis of the cluster data is based on the estimate of their 
r.m.s. peculiar velocity, U rms . The results are compared with predictions from cosmological models using 
linear theory. The statistics of the model velocity field is parameterized by the amplitude, % = ^s^if, 
and by the shape parameter, T, of the CDM-like power spectrum. This analysis is performed in redshift 
space, so as to circumvent the need to address corrections due to inhomogeneous Malmquist bias and 
to the redshift cutoff adopted in the sample selection. From the velocity correlation statistics we obtain 
7/s = 0.511q'o9 for r = 0.25 at the 2a level for one interesting fitting parameter. This result agrees with 
that obtained from a similar analysis of the SFI I-band Tully-Fisher (TF) survey of field Sc galaxies. 
Even though less constraining, a consistent result is obtained by comparing the measured U rms of clusters 
to linear theory predictions. For T = 0.25 we find r/g = O.BSIq' 22 , at la. Again, this result agrees, within 
the uncertainties, with that obtained from the SCI cluster sample based on TF distances. Overall, 
our results point toward a statistical concordance of the cosmic flows traced by spirals and early-type 
galaxies, with galaxy distances estimated using TF and D n -a distance indicators, respectively. 

Subject headings: Cosmology: observations - cosmology: theory - galaxies: distances and redshifts - 
large-scale structure of universe. 
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In this context, the recently completed all-sky redshift - 
distance survey of early-type galaxies (da Costa et al. 
2000a, hereafter ENEAR), probing a volume comparable 
to that of the existing catalogs of peculiar velocity data, 
is a welcome addition. The ENEAR galaxies sample dif- 
ferent regions of space and density regimes; the peculiar 
velocities are measured using an independent distance in- 
dicator; and the distances are based on separate types of 
observations, reduction techniques and corrections. Fi- 
nally, the ENEAR sample has well defined selection crite- 
ria, the completeness of the observations is uniform across 
the sky and the data, mostly new measurements by the 
same group, are in a homogeneous system. 

The present Letter has the twofold aim of comparing 
global statistical quantities, which describe the velocity 
fields traced by the TF and D n — a distance indicators, 
and of placing constraints on the nature of the fluctuation 
power-spectrum. Our analysis is based on the velocity cor- 
relation statistics and the r.m.s. one-dimensional peculiar 
velocity of clusters. These statistics were used by Borgani 
et al. (2000, BOO hereafter) and Borgani et al. (1997, B97 
hereafter) to analyze the SFI sample of field spirals and 
the SCI sample of cluster spirals (Giovanelli et al. 1997), 
respectively. In this paper, the same analysis is carried 
out for the ENEAR sample of field galaxies and groups 
and for ENEAR clusters (hereafter ENEARc; Bernardi et 
al. 2000, in preparation). 

2 THE DATA 

The ENEAR sample contains 1359 ellipticals brighter 
than rag = 14.5 with D n — a measured distances and 
569 cluster galaxies in 28 clusters (ENEARc). Galaxies 
have been objectively assigned to groups and clusters us- 
ing the information available from complete redshift sur- 
veys sampling the same volume. Our analysis is performed 
in redshift-space so as to avoid correcting for inhomoge- 
neous Malmquist bias and the redshift cutoff adopted in 
the sample selection. Therefore, we use the inverse D n — a 
template derived by Bernardi et al. (2000, in preparation) 
combining all the cluster data. Below, we limit our anal- 
ysis to objects within cz = 6000 km s -1 , so as to exclude 
those with very uncertain velocity measurements. This 
sub-sample consists of 355 field galaxies and 223 groups. 
In the cluster sample analysis we only consider the 20 clus- 
ters with cz < 6000 km s -1 . Of these, we discard the clus- 
ters CEN 30 and CEN 45; these systems lie along the 
same line-of-sight and are close in redshift-space making 
the assignment of galaxies to individual systems difficult. 
They are also suspected to form a bound system (Lucey & 
Carter 1988) and their large peculiar velocities observed 
may be due to non-linear effects. We also pay special at- 
tention to two other groups, AS714 and AS753, both with 
large peculiar velocities (~ 900 kms -1 ). These systems 
lie in the region of the Great Attractor and may also be 
subject to non-linear dynamical interactions. Below we 
discuss the impact of including or excluding these two sys- 
tems in the analysis. 



3 THE VELOCITY CORRELATION 
STATISTICS 

Our analysis of the velocity correlation statistics follows 
closely that presented in BOO. We refer to that paper for 
a more thorough discussion. We use the velocity correla- 
tion estimator originally introduced by Gorski et al. (1989, 
G89 hereafter): 
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where i?y is the angle between the direction of the i-th 
and the j-th object and the sums are over all the pairs at 
separation r in redshift space. In eq.(Q) Ui is the radial 
peculiar velocity of the i th object and we assign equal 
weight to all objects, so as to minimize the effect of cos- 
mic variance (see the discussion in BOO). The average 
of tpi(r) over an ensemble of cosmic flow realizations is 
*iW = (^iW) = A(r)"t>n(r) + [I- A(r)]V±(r), where 
'I'll and \&_|_ are the radial and transverse correlation func- 
tions of the three-dimensional peculiar velocity field (see 
G89). In linear theory, they are connected to the power- 
spectrum of density fluctuations, P(k), according to 
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where ji{x) is the i-th order spherical Bessel function and 
f(Q m ) ~ f^m 6 - The quantity A(r) is a moment of the 
selection function of the sample, which is fully specified 
by the spatial distribution of the objects in the sample 
(e.g., Gorski et al. 1989; BOO). Therefore, the model 
can be computed taking into account the specific sam- 
pling through the A(r) function. The velocity correla- 
tion function ipi(r) for the ENEAR sample is plotted in 
Figure 1 up to r — 3500 km s -1 , for all objects within 
cz = 6000 km s -1 . This separation range has been shown 
by BOO to be that where ipi is more stable for the SFI 
sample, which has about the same size as ENEAR. We 
choose the bin size to be 500kms _1 in order to keep these 
errors relatively small within each separation bin. We ver- 
ified that final constraints on the model parameter are left 
unchanged by halving the bin width. For the purpose of 
comparing ENEAR and SFI results, we show in Figure 
1 only the statistical errors due to the internal noise of 
the data set, which have been estimated as follows. At 
the position of each galaxy we add to the observed pe- 
culiar velocity a random component which is drawn from 
a Gaussian distribution having r.m.s. dispersion equal to 
the observational error reported for that object. Veloc- 
ity correlations are then computed for 1000 realizations of 
this perturbed data set and errors on ipi are estimated at 
each separation from the scatter among these realizations. 
Cosmic variance must not be included here, since ENEAR 
and SFI probe cosmic flows within the same region of the 
Universe. 
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Fig. 1. — The velocity correlation function, tpi(r) (in units of 
10 4 kms _1 ), for the ENEAR sample (filled circles). Open cir- 
cles and squares are the results for the SFI sample, as derived 
by BOO from two different zero-point calibrations of the TF re- 
lation. Errorbars, that for reasons of clarity are only reported 
for ENEAR, are the la uncertainties from the internal sample 
noise (see text). 

Remarkably the ipi velocity correlation of the ENEAR 
sample falls just between the two SFI estimates, based 
on the two zero-point calibrations of the TF relation pre- 
sented by BOO. This result contrasts with the disagreement 
originally found by G89 between spirals (Aaronson et al. 
1982) and ellipticals (Lynden-Bell et al. 1988). We use the 
ENEAR velocity correlation function to place constraints 
on cosmological models by following the same procedure 
discussed in BOO, and only briefly summarized here. We 
run N-body simulations for different cosmological models, 
and extract from each of them a fairly large number (256) 
of independent mock ENEAR samples. In each mock sam- 
ple, "galaxies" are placed in the same positions as in the 
real sample. The peculiar velocity of each mock galaxy 
is then perturbed with a Gaussian-distributed component 
associated to the observational error of its real counter- 
part. With this procedure, each set of mock samples in- 
cludes both cosmic variance and statistical noise. There- 
fore, we estimate the elements of the covariance matrix 

C t] = N~l ck j:lr ck - i>\) - where ^ 

is the value of the velocity correlation function at the i-th 
separation bin for the l-th mock sample, and ip\ is its en- 
semble average. Based on this approach, BOO showed that 

(a) linear theory provides a good description of the veloc- 
ity correlation statistics of the N-body simulated samples; 

(b) the relative amount of covariance, i.e. the values of 
C^/i/;\ipi, is independent of the underlying cosmology 

Based on these results, we compute a grid of linear- 
theory model predictions for ip\ (r) , as well as the elements 
of the corresponding covariance matrix expected for a sam- 
ple the size of ENEAR. We assume the power spectrum 
expression P(k) = AkT 2 (k), where the transfer function, 
T(k), is assumed to have the CDM-like form with the 
fc-dependence specified by the shape parameter T. The 
amplitude of P(k) is expressed in terms of <Tg, the r.m.s. 



fluctuation amplitude within a sphere of 8 /i _1 Mpc. There- 
fore, following eq.(Q), ipi (r) is entirely specified by the two 
parameters T and rjs — as^n- In order to derive con- 
straints on these parameters, we compute the weighted \ 2 
between the ENEAR correlation function, ipf NEAR , and 
that from model predictions, ^" lod , taking into account 
the covariance terms. The probability for model rejection 
is estimated, from the value of A\ 2 = X Xmim assuming 
a x 2 statistic, where Xmin ^ s the absolute minimum value. 

In Figure 2 we plot the iso-Ax 2 contours correspond- 
ing to lcr and 2cr confidence levels. The degeneracy of the 
constraint in the rjs-T plane is due to the fact that the 
coherence of the flow on a given scale depends not only 
on the overall amplitude of the power-spectrum but also 
on its slope. This is because peculiar velocities are gener- 
ated non-locally, so that coherence of the flow on a given 
scale can be associated either to fluctuations on compa- 
rable (large rjs and T) or on much larger scales (small rjg 
and r). Fixing T = 0.25, consistent with galaxy clus- 
tering data (e.g., Dodelson & Gaztanaga 1999), we find 
778 = 0.51+QQg at the 2<r level for one interesting fitting pa- 
rameter. We verified from the analysis of the mock samples 
that redshift -space distortions have a negligible effect in 
the estimate of tpi(r), with the resulting constraints on r/g 
being affected at most by about 5%. As expected from the 
comparison shown in Fig. 1, this result is in good agree- 
ment with that derived by BOO from the analysis of the SFI 
TF-survey of spiral galaxies. Therefore, we confirm that, 
for reasonable values of the power spectrum shape, the 
velocity correlation statistics favor small power-spectrum 
amplitudes. Although at variance with other analysis of 
velocity fields based on maximum likelihood analysis of 
the velocity correlation statistics (e.g. Zaroubi et al. 1998; 
Freudling et al. 1999; see the discussion in BOO), this re- 
sult agrees with independent constraints on the amplitude 
of the power spectrum, like those imposed by the num- 
ber density of local galaxy clusters (e.g., Eke et al. 1996; 
Girardi et al. 1998). 

4 The r.m.s. VELOCITY OF CLUSTERS 

The r.m.s. peculiar velocity of clusters has been used by 
several authors as further means to set constraints on cos- 
mological parameters (e.g., Moscardini et al. 1996; Bah- 
call & Oh 1996; B97; Watkins 1997). This analysis is 
repeated here for the ENEAR cluster sample. The obser- 
vational estimate from the ENEAR sample ranges from 
V r °£l = 450 ± 73 km s" 1 to V°£l = 470 ± 68 km s" 1 for the 
samples of 16 and 18 clusters, respectively, defined in Sec- 
tion || by either excluding or including AS714 and AS753. 
The error is the lcr scatter over 10 5 random realizations of 
the real sample, each generated from a Gaussian distribu- 
tion having the above T^ms and velocities convolved with 
the observational errors. From the theoretical side, linear 
theory for the growth of density fluctuations predicts that 
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the one-dimensional r.m.s. velocity is 
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where we use the expression W(kR) — exp(— k 2 R 2 /2) for 
the window function. Croft & Efstathiou (1994) verified 
that eq.(||) provides a rather good fit to results from N- 
body simulations for R values in the range 1.5-3 /i^ 1 Mpc. 
In the present analysis we adopt R — 1.5 h~ 1 Mpc, but 
point out that the results are largely insensitive to the ex- 
act choice. For instance, assuming R twice as large only 
increases the final constraints on rjs by about 8%. 




Fig. 2. — The la and 2a contours in the rjs-F plane from the 
analysis of the velocity correlation function, t[>i(r), for ENEAR 
groups and galaxies. Also shown with the dashed contours is 
the la confidence limits from the analysis of the r.m.s. peculiar 
velocity of ENEAR clusters. 

The procedure to establish the confidence level for a 
given model is the same as applied by B97. Let Vi and 
<7j be the velocity and its error for the i-th real cluster 
(i = 1, . . . , 16). For a given model, we generate Monte- 
Carlo samples, each containing 16 velocities, Vi, drawn 
from a Gaussian distribution, with dispersion provided by 
eq.(|^). For each sample, every cluster's velocity is esti- 
mated as a Gaussian deviate of the mean Vi and disper- 
sion Cj, and the r.m.s. velocity of the sample is then com- 
puted. For each model we generate N = 10 4 samples and 
then compute the fraction T with y/ ms (j = 1, . . . , N) at 
least as discrepant as V°^ with respect to their average 
value, Ylj ^-ms- Therefore, the smaller the value of 
T, the larger the probability for model rejection. After de- 
termining the highest value of T ', relative confidence levels 
are computed by determining standard decrements with 



respect to this maximum value (i.e., AJ 7 ~ 0.68 and 0.95 
for la and 2a exclusion levels). The resulting la con- 
straints on the r-?78 parameter space are shown in Fig. 2 
(dashed curves). Although this result is less constraining 
than that obtained from the velocity correlation analysis, 
it nicely overlaps with it, thus demonstrating that EN- 
EAR clusters and field galaxies consistently trace the same 
large-scale flows. For T = 0.25 we find rfs = 0.63 



+0.22 
-0.19 



at 



the la confidence level. This result is also consistent with 
that previously obtained from similar analyses of the SCI 
cluster velocities (B97, Watkins 1997). The inclusion of 
the AS714 and AS753 clusters in our analysis would only 
increase the resulting rj$ by about 5%. 

5 CONCLUSIONS 

We presented statistical analyses of the peculiar veloc- 
ity field within cz = 6000 km s -1 traced by field objects 
and clusters in the ENEAR sample based on D n -a dis- 
tances. We use the velocity correlation statistics tpi(r) 
to characterize the velocity field traced by field ellipti- 
cals and loose groups ad find results which are consistent 
with those obtained from the SFI sample of spirals with 
TF distances. Contrary to past claims, we find no statis- 
tically significant differences between the peculiar veloc- 
ity fields mapped by spirals and ellipticals. This result is 
in general agreement with and generalizes the findings of 
da Costa et al. (2000b) using the bulk-velocity statistics. 
Constraints on the power spectrum of density fluctuations 
were derived by resorting to linear theory. Assuming the 
shape of the power spectrum to be consistent with results 
from galaxy galaxy clustering analyses, T — 0.25, wc find 
77s = 0.5llg;Qg at 2a level for one interesting fitting pa- 
rameter. A consistent constraint is also obtained from 
the analysis of the r.m.s. velocity of ENEAR clusters; 
for the same value of the shape parameter T, it implies 
?/8 = 0.63^19 at thus consistent with results from 
the SCI cluster TF velocities (B97, Watkins 1997). Our 
results confirm the conclusion by BOO that the amplitude 
of cosmic flows can be reconciled with independent con- 
straints on the amplitude of density perturbations as that 
required by the number density of nearby rich clusters. 
They also show that consistent results are obtainable from 
independent distance indicators, once they are applied to 
homogeneously selected galaxy samples. 
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